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ABSTRACT 

Using density-functional theory with ultrasoft pseudopotentials, we previously 
investigated the structural and electronic properties of the low-pressure (cubic, tetragonal, 
and monoclinic) phases of Zr0 2 and HfC>2, in order to elucidate phonon modes, Born 
effective charge tensors, and especially the lattice dielectric response in these phases. We 
now extend this previous work by carrying out similar calculations on the two 
high-pressure orthorhombic phases, and by providing density-of-states and band-gap 
information on all polymorphs. Our results show that the electronic structures and 
dielectric responses are strongly phase-dependent. In particular, the monoclinic phases of 
Zt0 2 and Hf02 are found to have a strongly anisotropic dielectric tensor and a rather small 
orientational average (eb) compared to the two other low-pressure phases. Our calculations 
show that £q is even smaller in the orthorhombic phases. 



INTRODUCTION 

Dielectrics comprise a class of materials of great technical importance. The most 
widely- used dielectric in electronics is probably Si02, which has been persistently exploited 
in the computer industry for decades due to the excellent quality of the Si/Si02 interface. 
The dimensions of the gate-dielectric layer (currently Si02) have been rapidly reduced to 
the nanometer scale in order to enhance the speed of CMOS-related devices. Such 
reductions, however, have uncovered severe fundamental problems that imposed a physical 
limit on the use of Si02 films as gate dielectrics in the current Si/Si02 semiconductor 
technology. A great deal of effort has been directed to finding a way to overcome these 
limits. One particular direction is to use high-JT metal oxides with certain desired 
properties [1] to replace Si02 as the gate dielectric, since the higher dielectric constants of 
the oxides will allow one to make a physically thicker film while at the same time 
maintaining (or even increasing) the gate capacitance. 

Zr02, Hf02 and their structure-modified derivatives (e.g., Zr and Hf silicates) have 
emerged as important candidates for this purpose because they have much higher dielectric 
constants than Si02 (~20 compared with 3.5 for Si02) and because of their 
thermodynamic stability in contact with Si. Experiments have demonstrated their promise 
for this purpose (see Ref. [1] for a general review). It would be very desirable, therefore, to 
conduct a systematic theoretical study of the structural, electronic, and dielectric 
properties of Zr02 and Hf02, and of the dependence of these properties on the choice of 
crystalline polymorphs or possible modified structures. In two previous works, we studied 



the structural and vibrational properties of the three low-pressure phases (cubic, 
tetragonal, and monoclinic) of Zr0 2 [3] and Hf0 2 [4], especially the zone-centered phonon 
modes and the related dielectric properties. These previous works are extended here in two 
ways. First, we provide density-of-states (DOS) information on all polymorphs of Zr0 2 and 
Hf0 2 , which may be helpful in the interpretations of experiments on these materials. 
Second, we extend our previous calculations of dielectric properties to include the 
high-pressure orthorhombic phases using the ABINIT package [4]. We have also been 
engaged in simulating amorphous Zr0 2 (a-Zr0 2 ) using ab-initio molecular dynamics as 
implemented in the VASP [5] package, in a "melt-and-quench" fashion. The resulting 
amorphous Zr0 2 structures were analyzed in order to understand their local bonding 
geometry and other structural properties, and then used as inputs for a first-principles 
investigation of the vibrational and dielectric properties of the amorphous Zr0 2 systems 
using ABINIT. The results will be reported elsewhere. 

CRYSTAL STRUCTURES 

Zr0 2 (zirconia) is surprising similar to Hf0 2 (hafnia) in many physical and chemical 
properties. It is well established that both oxides exhibit three crystalline phases at 
ambient pressure, i.e., monoclinic, tetragonal and cubic, as temperature is increased. With 
increasing pressure, phase transitions involving two orthorhombic phases have been 
reported, although the exact structures of these high-pressure phases has not been 
established unambiguously ([6,7], and references therein). The first orthorhombic structure 
used in our simulation (conventionally labeled as Ortho-I) has space group Pbca; the other 
orthorhombic phase (Ortho-II) has space group Pnma. We study these phases at zero 
pressure, where they are apparently metastable. 

The cubic phase has the fluorite (CaF 2 ) structure, with 3 atoms (one chemical 
formula) per primitive unit cell. The tetragonal structure is a slight distortion of the cubic 
phase obtained by shifting alternating pairs of oxygen atoms up and down along z, and the 
number of atoms per primitive cell is doubled. The monoclinic phase, which has 12 atoms 
(four formulas units) per primitive cell, has a more complex structure and lower symmetry. 
In our calculations of the DOS functions, a 12-atom cell has been used for all of these 
phases, as in Refs. [2] and [3], in order to facilitate comparison. The monoclinic phase has 
two non-equivalent oxygen atoms, which evidently will have different site-projected DOS 
functions. Among the eight oxygens in the unit cell, four are threefold coordinated 
(denoted as Oi), while the other four (labeled as 2 ) are fourfold coordinated. Each 
threefold oxygen is bonded to three nearest-neighbor metal atoms in an almost planar 
configuration, while each four-fold oxygen forms a distorted tetrahedron with its four 
nearest metal neighbors. For the cubic and tetragonal phases, all oxygen atoms are 
equivalent, although we still use the notation Oi and 2 for convenience. Detailed 
discussion of the structures of these phases can be found in Refs. [2] and [3] . 

A primitive unit cell of the Ortho-I phase (for both Zr0 2 and Hf0 2 ) consists of 24 
atoms, while the Ortho-II phase has 12 atoms per unit cell. We used the lattice parameters 
given in [6] (Table III there) and [7] (Tables II and III there) as the starting point for our 



Table I: Structural parameters calculated for the orthorhombic phases of Zr0 2 and Hf0 2 . 
Lattice constants a, 6 and c, in A are given in the first row. Triplets of coordinates indicate 
the Wyckoff positions of the various atoms. 



Phase 


Zr0 2 


Hf0 2 


Ortho-I 


a=9.958, 6=5.224, c=5.006 
Zr (0.885, 0.035, 0.254) 

01 (0.792, 0.379, 0.133) 

2 (0.976, 0.739, 0.497) 


a=9.837, 6=5.129, c=4.954 
Hf (0.885, 0.035, 0.254) 

01 (0.791, 0.376, 0.128) 

2 (0.977, 0.739, 0.497) 


Ortho-II 


a=6.438, 6=5.521, c=3.273 
Zr (0.884, 0.256, 0.250) 

01 (0.164, 0.027, 0.250) 

2 (0.572, 0.142, 0.250) 


a=6.353, 6=5.435, c=3.240 
Hf (0.885, 0.256, 0.250) 

01 (0.163, 0.025, 0.250) 

2 (0.573, 0.141, 0.250) 



ground-state calculations for Ortho-I and Ortho-II, respectively. The Wyckoff coordinates 
define two non-equivalent sets of oxygen atoms, labeled by Oi and 2 . 

THEORETICAL METHODOLOGY 

The DOS functions are computed using a plane- wave implementation of density 
functional theory (DFT) in the local-density approximation (LDA) as parametrized by 
Ceperley and Alder [8], using ultra-soft pseudopotentials [9]. Although a 4x4x4 k-point 
sampling was found to provide sufficient precision in calculating total energies and forces, 
we used a refined mesh of k-points and higher 36- Ry cut-off to obtain smoother and more 
accurate DOS functions. The k-point meshes used in the DOS calculations are 20x20x20 
for cubic and tetragonal, 10 x 10x10 for monoclinic, 12x12x12 for Ortho-I, and 16x16x16 
for Ortho-II phases. In calculating the site-projected partial DOS, spherical radii of 0.9 A 
and 1.25 A are used for the metal ions (Zr or Hf) and oxygen ions respectively. For cubic 
and tetragonal phases, the larger 12-atom unit cell was used in the calculations, so the 
DOS values should be reduced accordingly if the DOS per primitive cell is desired. Further 
details can be found in [2] and [3]. 

The dielectric response of orthorhombic phases of Zr0 2 was calculated by specialized 
linear-response techniques as implemented in ABINIT [4], taking into account the coupling 
between phonons and the homogeneous electric field. A 4x4x4 k-point mesh, 
Troullier-Matins type pseudopotentials [10] (4 and 6 valence electrons for Zr and O 
respectively), and an energy cutoff of 35.0 Hartree provide satisfactory convergence. 

RESULTS 

In Table we present our calculated lattice parameters for the two orthorhombic 
phases, which agree well with the previous theoretical calculations [6,7]. Because the 
computed DOS functions for Hf0 2 and Zr0 2 are so similar, we only present those for Hf0 2 
here, shown in Figs. (l)-(3). Some comparisons are given later in Table HP for both oxides. 
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Figure 1: Density of states (DOS) for monoclinic HfC>2. Upper panels show the orbital- 
projected partial DOS for Hf, Oi, and O2, respectively. Lower panel shows the total DOS. 
Vertical dashed line indicates the valence-band maximum. 



Figure [U shows the total DOS (lowest panel) and the site-projected partial DOS (upper 
three panels) for monoclinic Hf02- We find this material to be an insulator with a band 
gap of 3.5 eV, in good agreement with previous theoretical results (3.48 eV in [11] and 
3.6 eV in [12]). The bands in the energy range from — 20 eV to — 15 eV are mostly due to O 
2s orbitals. The next bands, spanning between — 6eV and 0, are mainly O 2p in character, 
and Hf 5d orbitals dominate the unoccupied conduction bands ranging from 3 eV up to 
11 eV. The two non-equivalent oxygen sites have noticeably different DOS features, as 
expected from their different surroundings. 

In Fig. we provide the DOS functions calculated for cubic and tetragonal Hf02- It 
can be readily seen that these two phases have rather similar DOS features. Cubic Hf02 
has a band gap of 3.2 eV, while for tetragonal Hf02 it is 3.8 eV. Ref. [11] gives 3.4 eV and 
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Figure 2: Same as in Fig. 1, but now for cubic (a) and tetragonal (b) phases of Hf02- 
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Figure 3: Same as Fig. 1, but now for the two orthorhombic phases of HfC>2. 

3.77 eV for the cubic and tetragonal phases, respectively. Again, both theoretical works 
agree quite well. The most significant difference between these two phases and the 
monoclinic phase is that the Hf d-like states separate into two bands in the cubic and 
tetragonal phases, but form a single band in the monoclinic phase. The DOS functions for 
the two high-pressure orthorhombic phases are presented in Fig. from which one readily 
sees that the oxygen atoms exhibit two quite different types of DOS, as in the monoclinic 
phase. In fact, there is a surprising degree of resemblance of the monoclinic DOS with that 
of the Ortho-I and Ortho-II phases, especially the former. 

The DOS functions of the Zr02 polymorphs (not shown) are qualitatively very similar 
to the corresponding ones for Hf02- The most important differences are in the band gaps. 
Because of LDA error, the absolute values of the band gaps are not reliable, but trends 
should be meaningful. As shown in Table |Hj we find that the band gaps are systematically 
larger (by ~0.5 eV) for Hf02 than for Zr02, and that variations of crystal structure can 
lead to band gap changes of order 1 eV. Also tabulated in Table HU are our calculated 
dielectric constants (orientational average eb) for Zr02 and Hf02 polymorphs. For the 
low-pressure cubic, tetragonal, and monoclinic phases, the dielectric constants are quoted 
from [2] and [3] , where the full dielectric tensors are given. For the Ortho-I and II phases of 
Zr02, the linear-response scheme [4] is utilized in the calculation of the static dielectric 
tensors, which includes the lattice contributions as well as the electronic screening. In 
orthorhombic Zr02, eb becomes even smaller, approximately 20 and 18 for the Ortho-I and 
II phases, respectively. The dielectric tensors for the orthogrombic phases are diagonal, 
with elements (22.6, 18.1, 19.6) and (18.8, 18.9, and 17.8) for Ortho-I and Ortho-II 
respectively. The average principal values of the Born effective charge tensors are about 5.0 
and 2.5 for Zr and O atoms respectively for both orthorhombic phases. 

CONCLUSIONS 

We studied the structural, electronic, vibrational and dielectric properties for nearly 



Table II: Band gaps (E g ) and orientationally averaged dielectric constants (eo) calculated for 
crystalline phases of Zr0 2 and Hf0 2 . E g values in parentheses for Hf0 2 are from the theory 
of Ref. [11]. 



Phase 


ZrQ 2 


HfQ 2 




E g (eV) 


e"o 


E g (eV) 


eb 


Cubic 

Tetragonal 

Monoclinic 


2.63 
3.31 
2.98 


37 
38 
20 


3.15 (3.40) 
3.84 (3.77) 
3.45 (3.48) 


29 
70 
16-18 


Ortho-I 
Ortho-II 


3.06 
2.29 


20.1 
18.5 


3.75 
2.94 





all the currently well- recognized phases of Zr0 2 and Hf0 2 , using advanced first-principles 
techniques, with a special focus on the dielectric properties. It is found that the crystalline 
structure can have a substantial effect on the band gaps and dielectric constants. Our 
results show that the dielectric response depends strongly on the crystal phase, and can 
span a wide range of values (e = 18-40 for Zr0 2 ). Monoclinic Zr0 2 (or Hf0 2 ) has a 
rather anisotropic dielectric tensor and a smaller averaged dielectric constant than the 
tetragonal and cubic phase. The e~o of the Ortho-I and II phases appear even smaller. We 
hope that these theoretical results will shed light on the potential for Zr0 2 , Hf0 2 or their 
derivatives to replace Si0 2 as the gate dielectric in CMOS technology, and will provide 
hints for manipulating these oxides to serve this purpose. 
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